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Self/non-self discrimination is a fundamental requirement of life. Endogenous peptides presented by 
major histocompatibility complex class I (MHC I) molecules represent the essence of self for CDS T 
lymphocytes. These MHC I peptides (MIPs) are collectively referred to as the immunopeptidome. 
From a systems-level perspective, very little is known about the origin, composition and plasticity of 
the immunopeptidome. Here, we show that the immunopeptidome, and therefore the nature of 
the immune self, is plastic and moulded by cellular metabolic activity. By using a quantitative 
high-throughput mass spectrometry-based approach, we found that altering cellular metabolism via the 
inhibition of the mammalian target of rapamycin results in dynamic changes in the cell surface MIPs 
landscape. Moreover, we provide systems-level evidence that the immunopeptidome projects at the cell 
surface a representation of biochemical networks and metabolic events regulated at multiple levels 
inside the cell. Our findings open up new perspectives in systems immunology and predictive biology. 
Indeed, predicting variations in the immunopeptidome in response to cell-intrinsic and -extrinsic 
factors could be relevant to the rational design of immunotherapeutic interventions. 
Molecular Systems Biology 7: 533; published online 27 September 2011; doi:10.1038/msb.2011.68 
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Introduction 

While unicellular eukaryotes primarily employ self/non-self 
discrimination to avoid self-mating, multicellular organisms 
use self/non-self discrimination primarily in immune defense 
(Boehm, 2006} . Failure to respond to non-self can lead to death 
from infection whereas untoward response to self paves the 
way to autoimmunity. Peptides presented by major histocom- 
patibility complex class I (MHC 1} molecules represent the 
essence of self for CDS T lymphocytes (Rammensee etaU 1993; 
Klein et al, 2009). These MHC I-associated peptides (MIPs) 
regulate all key events that occur during the lifetime of CDS 
T cells in the thymus and the periphery (Goldrath and Bevan, 
1999; Klein et aU 2009). MIPs are presented at the surface 
of most nucleated cells in jawed vertebrates and are 
collectively referred to as the MHC I immunopeptidome 
(Istrail et aU 2004). 

Classic reductionist approaches have estabhshed the broad 
outlines of the MIP processing system and shown that 
MIPs derive from degradation of endogenous proteins by 
the proteasome and other peptidases (Yewdell et aU 2003; 
Hammer et al, 2007] . They have also highlighted the complex- 
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ity of the MIP repertoire. Thus, MHC I molecules present 
peptides encoded not only in the primary open reading frames 
but also those encoded in alternate reading frames (Starck and 
Shastri, 2011). Moreover, progresses in mass spectrometry 
(MS) have allowed increasingly sophisticated and compre- 
hensive large-scale analyses of MIPs (Mester et aZ, 2011). 
Large-scale analyses have yielded unprecedented insights into 
the peptide specificities and motifs of MHC molecules and the 
diversity of the MIP repertoire (Hunt et aU 1992; Engelhard 
et aU 1993; Bonner et aU 2002; Weinzierl et aU 200S}. They 
have also demonstrated that MIPs derive from all cell compart- 
ments and that the MIP repertoire can be modified by 
neoplastic transformation (Schirle et aZ, 2000; Hickman et al, 
2004; Weinzierl et aU 2007; Fortier et aU 200S}. Nevertheless, 
we still know very little about the genesis and molecular 
composition of the immunopeptidome: why do proteins such 
as STT3B yield abundant MIPs while others do not (Perreault, 
2010}? MS studies have also revealed that the immunopepti- 
dome is not a random sample of the proteome: many abundant 
proteins do not generate MIPs, while some low-abundance 
proteins generate large amounts of MIPs (Milner et aU 2006). 
Furthermore, large-scale analyses have yielded confUcting 
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results on the relation between the trans criptome and the MIP 
repertoire (Weinzierl et aU 2007; Fortier et aU 2008; Mester 
et aU 2011). Therefore, further systematic studies based on 
high-throughput technologies and integrative approaches are 
needed in order to elucidate the mechanisms that mould the 
immunopeptidome. In-depth mechanistic understanding of 
the immunopeptidome biogenesis would allow prediction of 
its molecular composition and would therefore be highly 
relevant to the development of immunotherapies (Zarling et aU 
2006; Sette and Rappuoh, 2010). 

A key unresolved question is whether the immunopepti- 
dome of a particular cell is plastic and affected by its metabohc 
state. This issue is of fundamental importance because cells 
targeted by CDS T lymphocytes (i.e. infected and neoplastic 
cells) are metabolically perturbed by intracellular parasite/ 
viruses or various transforming events. We addressed this 
question by studying the immunopeptidome of a mouse 
lymphoma cell line (EL4} treated with the classic mammalian 
target of rapamycin (mTOR) inhibitor, rapamycin. We selected 
this model because rapamycin is an exquisitely specific 
inhibitor of the serine/threonine kinase mTOR, and because 
mTOR is a central regulator of cellular homeostasis through its 
involvement in promoting anabolic and inhibiting catabolic 
processes (Sengupta etaU 2010} . Furthermore, mTOR signaling 
is regulated in numerous physiological and pathological 
conditions in all cell types, including neoplastic and immune 
cells (Araki et aU 2009; Caron et aU 2010; Sengupta et al 2010} . 

Results 

The effect of rapamycin on mTOR signaling 
in EL4 cells 

As a prelude to our peptidomic studies, we first assessed the 
effects of rapamycin on EL4 cells. The first mTOR complex, 
mTORCl, phosphorylates both S6K1 and 4E-BP1 (Dowling 
et aU 2010}. Recent studies have shown that rapamycin, 
previously thought to completely inhibit mTORCl activity, 
differentially affects 4E-BP1 and S6K1 (Choo et al 2008; Choo 
and Blenis, 2009; Thoreen et al 2009}. While rapamycin 
treatment completely and sustainably inhibits S6K1 activation, 
it partly and variably inhibits 4E-BP1 phosphorylation (Choo 
et al 2008}. To evaluate how rapamycin affected mTORCl 
signaling in EL4 cells, we treated EL4 cells with rapamycin for 
different time durations up to 48 h. As expected, Thr-389 on 
S6K1 was fully dephosphorylated upon rapamycin treatment, 
correlating with a decrease in cell size (Figure lA and B}. In 
contrast, Ser-65 and Thr-37/46 on 4E-BP1 were partially 
dephosphorylated by rapamycin (Figure lA}. As previously 
reported, we also noted that abundance of 4E-BP1 was 
decreased after 24 h of rapamycin treatment (Figure lA} 
(Dilling et al 2002}. Because 4E-BP1 is critically involved in 
cap-dependent translation via regulation of eIF4E, these 
results suggested that translation was maintained in EL4 cells 
in the presence of rapamycin. Accordingly, we observed that 
protein synthesis decreased during the first 12 h of rapamycin 
treatment, but progressively recovered thereafter (Figure IC}. 
Prolonged rapamycin treatment for up to 24 h has also been 
shown to inhibit the assembly of mT0RC2, the second mTOR 
complex (Sarbassov et al 2006}. Here, we observed that 
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mT0RC2 -mediated phosphorylation of ART at Ser-473 was 
transiently decreased but recovered and was increased 
after treatment for 48 h. Thus, rapamycin treatment for 48 h 
inhibited mTORCl, but activated mT0RC2 in EL4 cells. 
Collectively, these results showed that rapamycin-mediated 
mTORCl inhibition perturbed protein synthesis and cell size 
by differentially inhibiting S6K1 and 4E-BP1 in EL4 cells 
(Figure ID}. Of direct relevance to our analyses of the MIP 
landscape, rapamycin treatment did not affect expression 
levels of key proteins involved in the MHC I antigen processing 
and presentation pathway (Supplementary Figure SI}. 

High-throughput MS-based studies unveil the 
plasticity of the immunopeptidome 

The plasticity of the immunopeptidome can only be estimated 
by systems-level analyses (Benoist et al 2006; Germain et al 
2011}, and MS with high accuracy analyzers is the most 
comprehensive and versatile tool in large-scale proteomics 
(Yates et al 2009} . We therefore used a high-throughput MS- 
based approach (Fortier et al 2008; de Verteuil et al 2010} to 
profile in a time-sequential manner the abundance of 
MIPs presented by EL4 cells exposed to rapamycin for up to 
48 h (Figure 2A; Supplementary Information} . Subtraction of 
'contaminant peptides' eluted from (32m mutant EL4 cells 
allowed specific identification of genuine MIPs (Supplemen- 
tary Table SI} (Fortier et al 2008; de Verteuil et al 2010}. This 
resulted in identification and quantification of the relative 
abundance of 422 unique MIPs associated to different MHC I 
allehc products (Figure 2B; Supplementary Table S2}. The 
sahent finding was that the abundance of most peptides 
presented by classical MHC la molecules (H2D^ and H2K^} 
increased over time after 12-48 h of rapamycin treatment 
(Figure 2B}. In contrast, the abundance of most peptides 
associated to non-classical MHC lb molecules (Qal and Qa2} 
was not augmented after treatment with rapamycin for 48 h. 
Accordingly, exposure to rapamycin led to an increase in cell 
surface levels of H2D^ and H2K^ but not Qal and Qa2 
(Supplementary Figure S2B} . Clustering analysis of the time- 
series data sets revealed that the abundance of numerous MIPs 
increased progressively in response to rapamycin (Supple- 
mentary Figure S2A}. In general, the dynamic range of this 
response was modest with a median fold change in MIP 
abundance of 2.1 relative to untreated cells (calculated from 
Supplementary Table S2}. Nonetheless, a subset of 222 MIPs 
showed fold changes between 2 and 15, including 6 MIPs 
detected exclusively on rapamycin-treated cells (Supplemen- 
tary Table S2} . Perturbation of cellular metabolic activity with 
rapamycin therefore unveiled the plasticity of the immuno- 
peptidome, by increasing both the abundance and the diversity 
of MIPs at the cell surface in a time-dependent manner. 

mTOR inhibition induces functionally coherent 
changes in the transcriptome and the 
immunopeptidome 

Whether the immunopeptidome accurately projects intracel- 
lular metabolic changes to the cell surface is a fundamental 
question that has not been addressed yet. More precisely, we 
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Figure 1 Rapamycln differentially inhibits S6K1 versus 4E-BP1 in EL4 cells. Cells were treated with 20 ng/ml of rapamycln for the indicated time periods. (A) Levels of 
the indicated proteins were determined by western blotting. p-Actin served as a loading control. Data are representative of three independent experiments. (B) EL4 cells 
were treated or not (ctrl) with rapamycln for 48 h. One thousand cells were counted for each condition. Cell size was measured by light microscopy. The red line 
corresponds to the average cell size. Data are representative of three independent experiments. *P< 0.0005 (Student's Mest). (C) Relative protein synthesis was 
measured by pH]-leucine incorporation. Data (mean ± s.d.) are representative of three independent experiments. *P<0.01 , **P< 0.005 (Student's f-test). (D) Model for 
effects of rapamycin-mediated mTORCI inhibition in EL4 cells. Proteins and reactions were curated from the literature. Color code is based on the results in Figure 1 
after 48 h of rapamycln treatment. Activation and inhibition of components/reactions are depicted in green and red, respectively. Partial inactivation is represented by the 
thinner green lines. 



asked herein whether variations in the immunopeptidome 
reflected perturbations of rapamycin-sensitive metabohc 
processes regulated at the transcript level. To address this 
question, we first defined the impact of rapamycln on the 
transcriptome of EL4 cells. Using NimbleGen MM8 385K 
microarrays, we compared the transcriptome of untreated EL4 
cells to that of EL4 cells treated with rapamycln for 48 h. A total 
of 1353 differentially expressed transcripts were identified 
(fold change >2.0; P<0.05}, which represents about 3% of 
transcripts in EL4 cells (Figure 3 A) . Among these, 903 and 450 
mRNAs were under- and overexpressed, respectively. These 
genes are hereafter referred to as differentially expressed genes 
(DEGs). We next compared the 1353 DEGs with 98 unique 
source genes coding for the most differentially expressed MIPs 
(DEMs) (fold change >2.5; P<0.05: as described in Fortier 
et aU 2008} (Figure 3B; Supplementary Table S2}. Based on 
microarray gene expression data, only two DEM source genes 
(Tmodl and DhcrV) were found to be differentially expressed 
at the transcript level (Figure 3 A and B) . Interestingly, the only 
DEM that was less abundant was also found to be down- 
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regulated at the transcript level. Thus, these results indicate 
that only a minority of DEMs originate from DEGs. 

Analysis of Gene Ontology (GO) annotations revealed that 
101 and 70 genes (Figure 3C; Supplementary Figure S3B and C) 
coding for under- and overexpressed mRNAs were asso- 
ciated to 7 and 8 (Figure 3D} significantly enriched cellular 
processes, respectively (P<0.05} (Supplementary Table S3}. 
GO term analysis on the 98 DEM source genes also revealed a 
significant (P<0.05} enrichment for 6 cellular processes 
(Figure 3D} imphcating 38 DEM source genes (Figure 3C; 
Supplementary Figure S3 A; Supplementary Table S3}. Analy- 
sis of DEGs and DEM source genes involved in enriched 
cellular processes yielded several findings. First, there was no 
overlap between the 171 DEGs and the 38 DEM source genes 
(Figure 3C}. Second, cross-comparison of the overrepresented 
functional groups identified four distinct cellular processes 
(protein transport, cell cycle/proliferation, DNA rephcation 
and transcription} that were enriched in both DEGs and DEM 
source genes (Figure 3D}. In order to determine whether this 
result implies a significant functional relationship between 
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Figure 2 Rapamycin increases the abundance of MIPs presented by MHC la molecules. (A) Experimental design for identification and relative quantification of MIPs. 
EL4 cells were incubated with 20 ng/ml of rapamycin for 0 h (ctrl), 6, 12, 24 and 48 h. Cells from all experimental conditions were harvested simultaneously. MIPs from 
rapamycin-treated and Ctrl EL4 cells were isolated by mild acid elution and analyzed by nanoLC-MS-MS/MS. Heat maps displaying m/z, retention time and abundance 
were generated. A logarithmic intensity scale distinguishes between low (dark red) and high (bright yellow) abundance species. Analysis of p2m- mutant EL4 cells allows 
discrimination of genuine MIPs from contaminants (Fortier etal, 2008; de Verteuil etal, 2010). Examples of peptides that were differentially expressed (red line) or not 
(green line) between WT and p2m- mutant EL4 cells are highlighted in the boxes. (B) Volcano plots show the relative abundance of 416 H2^-associated peptides 
identified in three biological replicates at each time point. Six MIPs were detected uniquely after 48 h of rapamycin treatment. 



trans criptome and immunopeptidome variations, w^e devel- 
oped an all-pairs-shortest-path matrix, w^hich scores the 
functional connectivity betw^een tw^o hsts of genes 
(Figure 4A}. Using this matrix, w^e first calculated a con- 
nectivity score betw^een the 171 DEGs and the 38 DEM source 
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genes identified above. Then, bootstrapping w^as used as a 
statistical sampling method to calculate control connectivity 
scores from 10^ sets of 38 randomly selected MIP source genes 
from a database of 891 unique source genes encoding H2^- 
associated peptides (Supplementary Table S4} . These analyses 
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Figure 3 Rapamycin-mediated mTOR inhibition induces functionally coherent changes in the transcriptome and the immunopeptidome. (A) Volcano plot 
representation of the relative abundance of 42 586 transcripts after 48 h of rapamycin treatment. Transcripts were considered to be differentially expressed when the fold 
difference in abundance was >2.0 (P<0.05). Transcripts over- and underexpressed in rapamycin-treated cells relative to Ctrl cells are depicted in the red and green 
boxes, respectively. (B) Volcano plot representation of the relative abundance of 422 MIPs after 48 h of rapamycin treatment. MIPs that were the most differentially 
expressed (fold difference relative to untreated cells >2.5; P<0.05; fold difference and P-value based on Fortier et al, 2008) are highlighted in the blue boxes. 
(C, D) GO enrichment analyses were performed for DEM source genes and DEGs. In all, 38 DEM source genes (blue) were associated to 6 significantly enriched cellular 
processes. In all, 101 (green) and 70 (red) genes coding for under- and overexpressed mRNAs were associated to 7 and 8 significantly enriched cellular processes, 
respectively. (C) Venn diagram showing no overlap between DEM source genes and DEGs. Dashed boxes show representative genes that contributed to enrichment of 
four cellular processes in both DEM source genes and DEGs. (D) Venn diagram showing functional overlap between cellular processes overrepresented in DEM source 
genes and DEGs. The four cellular processes overrepresented in both DEM source genes and DEGs are listed in the dashed boxes. 



revealed that the 38 DEM source genes were tightly inter- 
connected to the 171 DEGs (bootstrapping; P=0.004} 
(Figure 4B} . Hence, while most DEMs did not originate directly 
from DEGs, they do originate from genes that are tightly 
functionally connected to DEGs. In other words, this systems- 
level analysis demonstrates that rapamycin-mediated mTOR 
inhibition induces functionally coherent changes in the 
transcriptome and the immunopeptidome. Given the tremen- 
dous complexity of the transcriptome, this finding suggests 
that the immunopeptidome of the cell is more complex than 
anticipated and that its plasticity might be extensive. 



DEMs arise from biochemical networks connected 
to mTOR 

The above results suggested that DEMs originated from genes 
connected to the mTOR network. However, considering that 
mTOR has a pervasive role in protein synthesis and degrada- 
tion (Caron et aU 2010), we could not discard the possibility 
that DEMs originated from some non-specific generic effect 
of rapamycin on protein metabohsm. In the latter case, 
DEM-coding genes would not be tightly connected to the 
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mTOR network. Hence, in order to further evaluate the 
relationship between DEMs and mTOR, we first conducted 
an analysis on the 98 DEM source genes (Figure 3B; 
Supplementary Table S2} using the interaction network 
database STITCH (Kuhn et al 2010). This analysis uncovered 
a network containing 30 DEM source genes that were 
interconnected and organized within discrete functional 
modules (Supplementary Figure S4}. Strikingly, the network 
included the chemicals rapamycin and everolimus (rapamycin 
analog) in addition to components (e.g. Rictor, Sgkl) and 
modules (e.g. mTOR signaling, translation, lipid biosynthesis) 
known to be directly regulated by mTOR (Caron et aU 2010). 
Therefore, we reasoned that rapamycin-mediated changes in 
the immunopeptidome might originate from genes that are 
very closely connected to components of the mTOR signaling 
network. To systematically evaluate this assumption, we 
measured the connectivity score between the 30 DEM source 
genes and components extracted from a comprehensive map 
of the mTOR interactome and signaling network (Supplemen- 
tary Table S5 based on Caron et aU 2010). By using the 
all-pair-shortest-path matrix described above, we calculated 
that the 30 DEM source genes were strongly interconnected 
to the mTOR network components relative to random 
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Figure 4 DEM source genes are tightly connected to transcriptonnic changes and the mTOR network. (A) An all-pairs-shortest-path matrix was developed by using 
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assignments (bootstrapping; P<10~^} (Figure 4C}. This 
systems-level analysis demonstrates that a substantial fraction 
of rapamycin-induced variations in the immunopeptidome 
arise from biochemical netw^orks w^hose components are 
highly connected to the target of rapamycin (i.e. mTOR). 
Thus, our data reinforce the notion that the immunopeptidome 
projects a functional representation of intracellular metabohc 
changes to the cell surface. Further immunopeptidomic 
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studies are needed to evaluate w^hether the mTOR netw^ork 
regulates the repertoire of MIPs in other contexts. 

To visualize the organization of relationships betw^een DEM 
source genes, the transcriptome and the mTOR netw^ork, 
w^e integrated in a global netw^ork the mTOR interactome and 
signaling netw^ork (Supplementary Table S5 based on Caron 
et aU 2010), the 98 DEM source genes (Figure 3B; Supplemen- 
tary Table S2} and the 171 DEGs (Figure 3C; Supplementary 
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Figure S3B and C). We found that 33 DEM source genes were 
constituents of 7 discrete functional subnetworks (Figure 5A} . 
All modules included components of the mTOR network. 
Integration of microarray data revealed that DEGs were found 



in 5 out of the 7 modules containing DEM source genes. 
Notably, two modules contained DEM source genes but no 
DEGs: the proteasome and the core mTOR signaling modules. 
Overall, 82% of the DEM source genes (27/33) were localized 
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Figure 5 DEM source genes are regulated at multiple layers within specific mTOR subnetworks. (A) The STITCH database (http://stitch.embl.de/) was used 
to generate a network from protein-protein interactions and functional associations involving identified DEM source genes, DEGs and mTOR network components. 
From the total network, subnetworks of DNA replication/transcription, cell cycle/proliferation, apoptosis, lipid biosynthesis, translation, proteasome and mTOR signaling 
were extracted. Legend for functional associations (edges) is depicted in Supplementary Figure S4. (B) EL4 cells were treated or not with rapamycin for 48 h. 
Total mRNA and NPR levels of the DEM source genes depicted in (A) were assessed by quantitative real-time PGR (see also Supplementary Table S6). 
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within subnetworks transcriptionally regulated by DEGs. 
Collectively, these data indicate that rapamycin-induced 
variations in the immunopeptidome arise mostly from 
transcriptionally perturbed subnetworks that are connected 
to components of the mTOR network. 



The immunopeptidome integrates events 
occurring at the transcriptional and 
CO- or post-translational level 

Generation of MIPs exploits generic cellular processes includ- 
ing transcription, translation and protein degradation (Yewdell 
etaU 2003; Fortier etaU 2008; de Verteuil etal 2010). However, 
the relative contribution of these different regulatory processes 
in moulding the immunopeptidome remains elusive. Here, we 
first sought to determine whether the 33 DEM source genes 
depicted in Figure 5A were translationally and/or transcrip- 
tionally regulated in response to rapamycin. After 48 h 
of rapamycin treatment, polysome-associated RNA as well as 
total RNA from both rapamycin-treated and -untreated cells 
were analyzed using quantitative real-time PGR. Translation of 
a large number of genes is reduced following rapamycin 
treatment (Grolleau et aU 2002). Accordingly, analysis 
of normalized polysomal mRNAs (NPR) (see Materials and 
methods) showed that, in rapamycin-treated cells, ribosomal 
loading efficiency was decreased for 36% (12/33), and 
increased for 6% (2/33) of the DEM source mRNAs 
(Figure 5B). In contrast, we observed that 36% (12/33) of 
the source genes were upregulated at the transcript level. 
Two points can be made from these data. First, while 
transcriptional activation could explain the increased abun- 
dance of about 36% of MIPs in rapamycin-treated cells, 
the role of translational activation appears negligible. As a 
corollary, increased levels of 20/33 MIPs can be ascribed to 
neither enhanced transcription nor translation of their source 
genes. From this, we infer that increased MIP levels on 
rapamycin-treated cells are caused mainly by co- and/or post- 
translational mechanisms. 

In order to probe the role of co- and/or post-translational 
mechanisms, we first assessed the abundance and stability of 
eight DEM source proteins by western blotting on total lysates 
from EL4 cells treated or not with rapamycin and cyclohex- 
imide (Figure 6; Supplementary Table S6). The abundance of 
two proteins was increased (Myb and Tfdp2) while that of the 
other six proteins was unchanged (Rfc2, Psmb9, Psmb8, 
Pabpcl, Rictor and Sgkl). Notably, the eight DEM source 
proteins were relatively stable with a half-hfe ^8h in the 
presence or absence of rapamycin. Myb was the only DEM 
source protein degraded more rapidly (two-fold; half-hfe from 
8 to 4h) in the presence of rapamycin. These data strongly 
suggest that MIP levels do not correlate with the abundance or 
stability of the native form of DEM source proteins. However, 
up to 30 % of proteins are so rapidly degraded (with a mean ti/2 
of ^lOmin) by the ubiquitin-proteasome that they escape 
detection in standard cycloheximide chase experiments 
(Qian et aU 2006). Most of these rapidly degraded proteins 
are defective ribosomal products (DRiPs) resulting from 
imperfections of protein synthesis or assembly (Yewdell and 
Nicchitta, 2006) . Accordingly, most protein species have two 
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half-lives: a longer one for the well-conformed native protein 
and a short one for DRiPs. Importantly, MIPs are generated at 
higher efficiency from rapidly degraded proteins than from old 
proteins (Yewdell et al, 2003; Yewdell and Nicchitta, 2006). 
Western blotting experiments on whole EL4 cell extracts 
revealed that rapamycin increased the amount of proteasomal 
substrates, that is, ubiquitinated proteins that accumulate in 
lactacystin-treated cells (Figure 7A). We focused on Rictor, a 
DEM source protein that is a component of the mTOR complex 
2, to evaluate more specifically whether rapamycin enhances 
proteasomal substrates generation from DEM source proteins. 
Rictor is coding for the KALSYASL peptide, which was 
exclusively detected after 48 h of rapamycin treatment 
(Supplementary Table S2) . EL4 cells were treated or not with 
rapamycin for 48 h, with or without lactacystin for the last 8 h. 
In the presence of lactacystin, we observed a substantial 
accumulation of polyubiquitinated Rictor proteins in rapamy- 
cin-treated cells relative to untreated cells (Figure 7B). This 
result indicates that more polyubiquitinated Rictor proteins 
are subjected to proteasome-mediated degradation upon 
rapamycin treatment. Because steady-state levels of the native 
Rictor protein were not affected by rapamycin (Figure 6) , we 
infer that the overabundant Rictor MIP derives from a rapidly 
degraded pool of Rictor and not from the long-lived native 
form of Rictor. Gollectively, our data indicate that changes in 
the immunopeptidome integrate events occurring at the 
transcriptional and co- or post-translational level. Further- 
more, consistent with the fact that mTOR links protein quality 
and quantity control (Qian et aU 2010), our work suggests that 
enhanced generation of rapidly degraded proteasomal sub- 
strates is instrumental to rapamycin-induced changes in the 
immunopeptidome. Altogether, our results show that the 
immunopeptidome projects at the cell surface a representation 
of biochemical networks and metabohc events regulated 
at multiple levels inside the cell. 



MIPs appearing de novo on rapamycin-treated 
cells elicit cytotoxic T-cell responses 

Finally, we wished to determine whether MIPs that were 
exclusively detected on rapamycin-treated cells could be 
antigenic. To test this, we selected the KALSYASL and the 
VNTHFSHL peptide, which are encoded by Rictor and Ttc21b, 
respectively (Supplementary Table S2). We immunized 
G57BL/6 mice with dendritic cells (DGs) coated with 
KALSYASL or VNTHFSHL synthetic peptide. Splenocytes from 
immunized mice did not kill untreated EL4 cells, but they 
showed specific cytotoxicity for EL4 cells treated with 
rapamycin or EL4 cells coated with the peptide used for T- 
cell priming (Figure 7G and D; Supplementary Figure S5). 
Further studies are needed to evaluate (1) the proportion of 
novel or overabundant MIPs that are antigenic and (2) whether 
particular pathways and cellular processes generate MIPs that 
are more likely to be antigenic. Nevertheless, the antigenicity 
of KALSYASL and VNTHFSHL is consistent with the notion 
that the immune system is tolerant to MIPs expressed 
at physiological levels but can mount biologically relevant 
immune responses toward self MIPs present in excessive 
amounts (Schild et aU 1990). 
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Figure 6 Relative abundance and stability of eight DEM source proteins in the presence or absence of rapamycin. EL4 cells were treated or not (ctrl) with 20 ng/ml of 
rapamycin for 48 h. Levels of DEM source proteins were determined by western blotting and quantified by densitometry. Fold change indicates the rapamycin/ctrl protein 
abundance ratio. To estimate the half-life of DEM source proteins, cells were treated with cycloheximide (CHX) for 0, 4 and 8 h in the presence and absence of 
rapamycin. Results in the graphs are expressed as a percentage of the remaining protein abundance in CHX-treated cells relative to CHX-untreated cells. p-Actin served 
as a loading control. One representative western blot out of three is shown for individual proteins (P-value: Student's f-test). Source data is available for this figure in the 
Supplementry Information. 
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Figure 7 Rapamycin-treated cells contain increased levels of proteasomal substrates and express antigenic MIPs. (A, B) EL4 cells were treated or not with 20 ng/ml 
rapamycin for 48 h, combined or not with 10 |iM of lactacystln for the last 8 h. (A) Accumulation of ubiquitinated proteins was measured on total cell lysates by western 
blotting. (B) Polyubiquitinated proteins from total cell lysates were isolated with TUBEs and levels of polyubiquitinated Rictor were then determined by western blotting. 
One representative western blot out of three is shown. (C, D) Mice were immunized with DCs coated with VNTHFSHL (C) or KALSYASL (D) peptide. Splenocytes from 
primed mice were tested for cytotoxic activity against CFSE-labeled target EL4 cells at different E/T ratios. EL4 cells coated with VNTHFSHL or KALSYASL were used as 
positive control. Data represent the mean ± s.d. for three mice per group. Source data is available for this figure in the Supplementry Information. 



Discussion 

Early proteomic studies were conducted with analyzers whose 
sensitivity (dynamic range) and accuracy were orders of 
magnitude inferior to that of MS analyzers that are now 
available (Depontieu etaU 2009; Yates etaU 2009; Nilsson etaU 
2010). As a result, early studies on the immunopeptidome 
identified only the more abundant MIPs. Those MIPs were 
found to derive from highly abundant housekeeping proteins 
that are common to many cell types (Marrack et aU 1993; 
Hughes and Hughes, 1995; Barnea etaU 2002; Engelhard etaU 
2002). More recently, high-throughput MS-based analyses 
have shown that the immunopeptidome conceals a cell-type- 
specific signature. Thus, though the immunopeptidome of DCs 
and thymocytes partially overlap, no less than 40% of their 
MIPs are cell-type specific (Fortier et aU 2008; de Verteuil et aU 
2010} . The present study demonstrates that perturbation of a 
single signaling pathway can lead to significant changes in the 
composition of the immunopeptidome. More specifically, 
our work shows that dynamic changes of mTOR signals are 
reflected in MHC I presentation of numerous peptides 
associated with the mTOR interactome and its signaling 
network. An important question is whether changes in the 
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MIP repertoire induced by rapamycin would be found in 
response to any type of cell stress. Two elements strongly 
support our contention that rapamycin-induced changes in the 
immunopeptidome are connected to the mTOR network. First, 
we found that the connectivity between DEM source genes and 
mTOR network components was amazingly strong (P< 10~^}. 
A bootstrap procedure (500 000 iterations) failed to reveal a 
single set of peptide source genes that were so tightly 
connected to component of the mTOR network (Figure 4C}. 
Second, we have previously evaluated the impact of other 
types of cell stress (tunicamycin, palmitate or glucose 
deprivation) on expression of MHC-peptide complexes at the 
surface EL4 cells (the cell line used the present study) 
(Granados et aU 2009). We found that cell stress induced 
by tunicamycin, palmitate and glucose deprivation led to a 
decrease in the generation of MIPs and in the expression of 
MHC I molecule. These effects are the exact opposite of what 
we found upon treatment of EL4 cells with rapamycin: 
overexpression of MHC I molecules (Supplementary Figure 
S2B) and increased abundance of cell surface MIPs 
(Figure 2B). In our opinion, it is therefore sound to conclude 
that rapamycin-induced variations in the immunopeptidome 
(1) arise from biochemical networks whose components are 
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highly connected to the target of rapamycin (i.e. mTOR) and 
(2) are not unspecific changes found with any type of cell 
stress. However, mTOR senses and integrates multiple 
environmental cues and mTOR network components interact 
with many signaling pathways, such as the Wnt pathway 
(Inoki et aU 2006; Caron et aU 2010} . Therefore, further studies 
will be needed to determine whether perturbation of mTOR by 
other agents or in other contexts would bring the same changes 
in the MIP repertoire as those induced by rapamycin. 

An important imphcation for the immune system is that, at 
the peptide level, the immune self is plastic and its molecular 
composition is influenced by the cell's metabohc activity. We 
found that variations in the immunopeptidome functionally 
reflected perturbations of rapamycin-sensitive metabohc 
processes regulated in the transcriptome. More precisely, we 
observed that changes in the MIP repertoire originate from 
genes that are very closely functionally connected to DEGs. 
We observed that only 2 % of DEMs did originate directly from 
DEGs as measured by DNA microarray experiments. In 
contrast, quantitative PGR revealed that 36% of the DEM 
source genes were slightly upregulated at the transcript level. 
In fact, eight DEM source genes (Bptf, Bub lb, Ccnf, Chaflb, 
Dtl Mnath Myb and TfdpZ) in Figure 3 come from DEGs that 
were identified by quantitative PGR (fold change between 1.2 
and 3.1; Supplementary Table S6} but missed by microarrays. 
Differences in the sensitivity of both techniques are likely to 
explain this discrepancy. Indeed, whereas microarrays permit 
genome-wide profiling of mRNA expression levels, traditional 
quantitative PGR is by far more sensitive (Germain et aU 2011} . 
Therefore, genome-wide quantitative PGR would have 
revealed a more significant overlap between DEM source 
genes and DEGs. Although the overlap between DEM source 
genes and DEGs was underestimated from microarray data, 
altogether our data indicate that changes in MIP abundance 
originate mostly from post-transcriptional mechanisms. 
Furthermore, data integration analysis revealed that 82% of 
DEMs arised from subnetworks that were transcriptionally 
reorganized in response to rapamycin. Thus, modulation in 
the composition of the immunopeptidome might originate 
predominantly from gene-regulatory subnetworks that are 
transcriptionally perturbed in response to multiple cellular 
stresses (i.e. pathogen infection (Kumar et aU 2010}}. In other 
words, genes involved in transcriptionally perturbed subnet- 
works, although not necessarily differentially expressed at the 
transcript level, might participate in modulating the composi- 
tion of the immunopeptidome. Gonsequently, identification of 
components within transcriptionally perturbed subnetworks 
might enable prediction of peptide source genes modulating 
the immunopeptidome's composition. 

mTOR integrates environmental cues in terms of nutrients 
and growth factors. Therefore, it is reasonable to infer that cell- 
extrinsic factors modulate the repertoire of self peptides 
presented by MHG I. In addition, our results suggest that the 
peptide MHG I repertoire is a functional, if not a direct, 
representation of the transcriptional state of a cell. Together, 
that means the nature of the immune self is much more 
complex than anticipated. Moreover, the notion that MIPs 
presented specifically on metabohcally stressed cells can be 
antigenic (e.g. KALSYASL and VNTHFSHL} could have several 
imphcations that need to be explored: immune responses 
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against metabohc stress-associated MIPs could contribute to 
elimination of infected or transformed cells (Gleimer and 
Parham, 2003}, but might also ehcit autoimmunity (Todd etaU 
2008} . In keeping with this hypothesis, some self MIPs were 
found to be upregulated on infected cells, and some of them 
were immunogenic (Herberts etaU 2003; Hickman etaU 2003; 
Ovsyannikova et aU 2005} . 

From a more general perspective, our data provide a proof of 
concept that acquisition and integration of large-scale, 
quantitative biological data from multiple regulatory layers 
(i.e. -omics} may enable prediction of variations in the 
composition of the immunopeptidome. The notion that the 
immunopeptidome is plastic and conveys at the cell surface an 
integrative view of cellular protein metabohsm opens up new 
perspectives in systems immunology and predictive biology 
(Benoist et aU 2006; Germain et aU 2011} . The immune system 
has co-evolved with pathogens and positive selection imposed 
by pathogens has hastened the evolution of immune genes and 
thereby increased the complexity of the immune system 
(Waterston et aU 2002; Hedrick, 2004}. Given the complexity 
of the immune system and the quintessential importance of 
self/non-self discrimination, it is imperative to further develop 
and exploit systems-level quantitative methods that will 
enable modeling of the immunopeptidome's plasticity. Such 
systems-oriented approaches will help to simulate and predict 
how cell-autonomous and -extrinsic elements modulate the 
composition of the immunopeptidome in health and disease. 

Materials and methods 

Cell culture, flow cytometry and determination 
of cell size 

EL4 cell lines were cultured as previously reported (Fortier et al, 2008} . 
Rapamycin, lactacystin and cycloheximide were purchased from 
Sigma-Aldrich. For MHC I labeling, cells were stained with anti- 
H2D^ (B22-249.R1; Cedarlane}, anti-H2K^ (Y3; ATCC}, anti-Qa2 
(clone 1-1-2; BD Pharmingen), anti-Qalb (6A8.6F10.1A6) antibodies 
and analyzed on a BD LSR II flow cytometer using FACSDiva software 
(BD Biosciences) (Meunier et al, 2005) . Cell number and cell size were 
measured by light microscopy using an Axio Imager microscope 
(Zeiss) and the Northern Eclipse software (Empix Imaging Inc.). 



Western blotting and isolation of 
polyubiquitinated proteins 

EL4 cells were harvested and lysed in RIPA buffer, samples were 
resolved by SDS-PAGE and immunoblotted with the following 
antibodies: ERAAP (gift from N Shastri, University of California, 
Berkeley, CA); Tapasin (gift from TH Hansen, Washington University 
School of Medicine, St Louis, MO); Lmp7, Lmp2, pi, pS, ERp57, Sgkl, 
P2m, PDI (Abeam); Calnexin, Calreticulin, p-actin (Sigma); PA28, 
S6K1, P-S6K1 (Thr-389), 4E-BP1, p-4E-BPl (Ser-65), p-4E-BPl (Thr- 
37/46), a-tubulin, Rictor, Akt, p-Akt (S473), Ubiquitin (Cell Signaling); 
Myb (Upstate, Millipore); MecU, P2 (Biomol International); Rfc2 
(Proteintech Group, PTG); Tfdp2 (Proteintech Group, PTG); Pabpcl 
(Cell Signaling); H2KVH2D^ (2G5; Thermo Scientific); Tapl, Tap2 
(Santa Cruz Biotechnology). Isolation of polyubiquitinated proteins 
was achieved by using Tandem Ubiquitin Binding Entities (TUBEs) 
(LifeSensors, Inc.) according to the manufacturer's instruction. 



Metabolic radiolabeling 

To measure protein synthesis, EL4 cells were cultured in the absence or 
presence of rapamycin (20ng/ml) for Oh (ctrl), 6, 12, 24 and 48 h. 
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During the last hour, cells were transferred to 96- well plates (0.5 x 
10^ cells/well, 12 replicates per conditions) and [^H]-leucine (10|iCi/ml) 
(Perkin-Elmer) was added to each well. At the end of the culture period, 
cells were harvested with a Cell Harvester (Perkin-Elmer) on a fiber glass 
filter (Packard) for washing. Filters were dried for 24 h and 20|il of 
scintillation liquid per well (MicroScint-0, Perkin-Elmer) was added. 
Radioactivity incorporation was determined with the TopCount-NXT 
microplate scintillation and luminescence counter (Packard) . 



MIP classification and abundance profiling 

Each MIP listed in Supplementary Tables SI and S2 was classified 
according to restriction size and binding motif (Fortier et al, 2008; de 
Verteuil et al, 2010) . The raw ratios for the time profiles of detected 
MIPs were then clustered using the TIGR MultiExperiment Viewer 
(TMEV) software. Hierarchical clustering with Euclidean distances 
and average linkage was used. 



Peptide extraction and MS analyses 

EL4 cells were treated with 20 ng/ml of rapamycin for 0, 6, 12, 24 and 
48 h. MIPs were then isolated from EL4 cells as described previously 
using mild acid elution (Fortier et al, 2008; de Verteuil et al, 2010). 
Three biological replicates (5 x 10^ EL4 cells per replicate) were 
prepared. MIPs obtained after acid elution were separated using an 
off-line 1100 series binary LC system (Agilent Technologies) to remove 
contaminating species. Peptides were loaded on a homemade strong 
cation exchange (SCX) column (0.3 mm internal diameter x 50 mm 
length) packed with SCX bulk material (Polysulfoethyl ATM, PolyLC). 
Peptides were fractionated with a gradient of 0-25% B after 25min, 
25-60% B after 35min (solvent A=5mmol/l ammonium formate, 
15% acetonitrile, pH 3.0; solvent B=2mol/l ammonium formate, 
15% acetonitrile, pH 3.0). MIPs were collected in five consecutive 
fractions and brought to dryness using a speedvac. MIP fractions were 
resuspended in 2% aqueous acetonitrile (0.2% formic acid) and 
analyzed by nanoLC-MS/MS on a LTQ-Orbitrap mass spectrometer 
(Thermo Fisher Scientific) (Fortier et al, 2008). Full mass spectra were 
acquired with the Orbitrap analyzer operated at a resolving power of 
60 000 (at m/z 400) and colhsion-activated dissociation tandem mass 
spectra were acquired in data-dependent mode with the quadrupole 
linear ion trap analyzer. Mass calibration used an internal lock mass 
(protonated (Si(CH3)20))6; m/z 445.12057) and typically provided 
mass accuracy within 5p.p.m. for all nanoLC-MS/MS experiments. 

MS/MS sequencing, hierarchical clustering and 
identification of MIP source proteins 

Peptidomic data were analyzed using Xcalibur software and peak lists 
were generated using Mascot distiller (version 2.1.1, http://www. 
matrixscience.com). Database searches were performed against the 
International Protein Index mouse database (version 3.23 containing 
51 536 sequences and 24 497 860 residues) using Mascot (version 2.2, 
http://www.matrixscience.com) with a mass precursor tolerance of 
±0.05 Da and a fragment tolerance of ±0.5 Da. Searches were 
performed without enzyme specificity and a variable modification of 
oxidized Met. All search results were filtered using an MHC motif filter 
based on the predicted mouse MHC I allele motifs (Supplementary 
Information). Raw data files were converted to peptide maps 
comprising m/z values, charge state, retention time and intensity for 
all detected ions above a threshold of 10 000 counts using in-house 
software (Mass Sense) (Fortier et al, 2008; de Verteuil et al, 2010). 
Peptide maps were aligned and clustered together to profile the 
abundance of Mascot identified peptides using hierarchical clustering 
with criteria based on m/z and time tolerance (±0.02 m/z and 
± 1 min) . This resulted in a list of non-redundant peptide clusters for 
all replicates of all samples to be compared. Peptides found in eluates 
from P2m mutant EL4 cells were considered 'contaminant peptides'. 
Subtraction of contaminant peptides eluted from untreated or 
rapamycin treated (for 48 h) P2m EL4 cells allowed specific identifica- 
tion of genuine MIPs (Fortier et al, 2008; de Verteuil et al, 2010). MIPs 
were further inspected for mass accuracy and MS/MS spectra were 
validated manually. The Sidekick resource (http://www.bioinfo. 
iric.ca/sidekick/Main) was used to identify MIP source proteins. The 
list of MIPs reported in the present work has been provided to The 
Immune Epitope Database and Analysis Resource (http://beta.immu- 
neepitope.org/) (Vita et al, 2010). MS/MS of all peptide identifications 
are available at http://www.thibault.iric.ca/proteoconnections 
(Courcelles et al, 2011) under the project name: MHC I immunopepti- 
dome, and were submitted to the PRIDE database (Vizcamo et al, 2010) 
under the accession numbers: 18855, 18856. 



Microarray analysis 

EL4 cells were treated or not with 20 ng/ml of rapamycin for 48 h. 
Three biological replicates were prepared. Total RNA was extracted 
from EL4 cells with TRIzol RNA reagent (Invitrogen) as instructed by 
the manufacturer. Samples were purified using DNase (Qiagen) and 
the RNeasy Mini kit (Qiagen), and the overall quality was analyzed 
with the 2100 Bioanalyzer (Agilent Technologies). Purified RNA 
(10 |ig/sample) was hybridized on MM8 385K NimbleGen chips 
according to the manufacturer's instruction. Arrays were scanned 
using a GenePix4000B scanner (Axon Instruments) at 5 \im resolution. 
Data were extracted and normalized using the NimbleScan 2.4 
extraction software (NimbleGen Systems). Further microarray ana- 
lyses were performed using GeneSpring GX 7.3.1. The complete 
microarray data sets have been deposited in ArrayExpress (http:// 
www.ebi.ac.uk/arrayexpress) under accession number E-MEXP-3007. 

Isolation of total and polysomal RNA for 
quantitative real-time PCR 

EL4 cells were treated or not with 20 ng/ml of rapamycin for 48 h. 
Three biological replicates were prepared. Total RNA was extracted as 
described above. Polysomal RNA was extracted from 2 x 10^ EL4 cells 
as described previously (Rajasekhar et al, 2003). Gene expression 
level was determined using primer and probe sets from Universal 
ProbeLibrary (https://www.roche-applied-science.com/sis/rtpcr/ 
upl/index.jsp). Primer sequences are given in Supplementary Table 
S7. PCR reactions for 3 84- well plate format were performed as 
described in Baron et al (2007) . All reactions were run in triplicate, and 
the mean values were used for quantification. The mouse Hprtl and 
Tbp genes were used as endogenous controls for total RNA and 
polysomal RNA, respectively. In order to obtain an unbiased profile of 
the normalized ribosomal loading efficiencies for cells, we calculated 
the net polysomal mRNA to total mRNA ratios. Such data are presented 
in the manuscript as NPR (Rajasekhar et al, 2003). 



GO enrichment and connectivity score 

The Database for Annotation, Visualization and Integrated Discovery 
(DAVID) Bioinformatics Resource (Peters and Sette, 2007) was used to 
identify significantly enriched GO terms associated to DEM source 
genes and DEGs from untreated and rapamycin-treated EL4 cells. 
To calculate connectivity scores between two lists of nodes (genes/ 
proteins), computed scores from STRING (http://string-db.org/) 
(Szklarczyk et al, 2011) were used. By using computed scores in 
STRING, functional associations between two lists of nodes were 
integrated into an all-pairs-shortest-path distance matrix. Each 
functional association in the all-pairs-shortest-path matrix was 
transformed into a 'distance' (D), defined as -ln(S), where S is the 
computed STRING score. The all-pairs-shortest-path matrix contains 
the length of the shortest path (distance) between every pair of nodes 
in the network. A connectivity score was then obtained by calculating 
the mean of the shortest path distance between every pair of nodes in 
a given matrix. 



STITCH network 

The 98 DEM source genes highlighted in Figure 3B were used to build 
the DEM network (Supplementary Figure S4). The 98 source genes 
encoding DEMs, the 171 DEGs classified within the enriched GO terms 
(P<0.05) in Supplementary Figure S3B and C and the mTOR network 
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components in Supplementary Table S5 were used to build the global 
rapamycin-regulatory network. Subnetworks (Figure 5A) were manu- 
ally extracted based on the modular organization of the total network 
(Hartwell et al, 1999; Rives and Galitski, 2003). Each of these 
molecules was individually searched in the STITCH (http://stitch. 
embl.de/) (Kuhn et al, 2010) database for its functional associations. 
The parameters used to elucidate the network were as follows: 
40 additional nodes, a network depth=l, and interactions with a 
minimum STITCH combined score of 0.400. Functional associations 
identified for all of the molecules were then imported into Adobe 
Illustrator as a Scalable Vector Graphics image (action view 
configuration) . 



In vitro cytotoxicity assay 

Bone marrow-derived DCs were generated as previously described 
(Fortier et al, 2008; de Verteuil et al, 2010) . On day 9 of culture, peptides 
were added at a concentration of 2 x 10~^ M and incubated with DCs 
for 3h at 37 °C. For mouse immunization, 10^ peptide-pulsed DCs 
were injected i.v. in C57BL/6 females at day 0 and day 7. On day 14, 
splenocytes were harvested from the spleens of immunized mice and 
depleted of red blood cells using 0.83% NH4CI. Cells were plated 
at 5 X 10^ cells/well in 24-well plates and restimulated with 2 x 10~*^M 
peptide at 37 °C. After 6 days, cytotoxicity was evaluated by a 4-h 
CFSE-based assay. The percentage of specific lysis was calculated as 
follows: (number of remaining CFSE + cells after incubation of target 
cells alone - number of remaining CFSE + cells after incubation with 
effector cells) /number of CFSE+ cells after incubation of target cells 
alone x 100. 



Supplementary information 

Supplementary information is available at the Molecular Systems 
Biology website (www.nature.com/msb). 
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